Calorons and localization of quark eigenvectors in lattice QCD 
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We analyze the localization properties for eigenvectors of the Dirac operator in quenched lattice 
QCD in the vicinity of the deconfinement phase transition. Studying the characteristic differences 
between the Z3 sectors above the critical temperature Tc, we find indications for the presence of 
calorons. 
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One of the most discussed topics in hadron physics is 
the chiral phase transition of QCD, and the microscopic 
processes connected with it. Many current and planned 
experiments are at least partially motivated by the hope 
that they will shed some light on the issues involved. 

In this contribution we study the localization proper- 
ties of the lattice Dirac operator in the neighborhood of 
the chiral phase transition. We are motivated by the 
fact that one of the most popular pictures of the phase 
transition relates it to the properties of QCD instan- 
tons and anti-instantons (see e.g. Ref. |l|). The origin 
of this connection is the following observation: For each 
isolated instanton or anti-instanton there exists a local- 
ized zero mode of the Dirac operator. For a liquid of 
instantons and anti-instantons these zero modes should 
be perturbed to form a band of small eigenvalues. At 
higher temperatures it is thought that instantons and 
anti-instantons may pair to form molecules, and that the 
associated modes will no longer have particularly small 
eigenvalues, but instead become an inconspicuous part of 
the bulk spectrum. 

The Banks-Casher formula 101, 
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relates the chiral condensate {qq) to /3Dirac(0), the den- 
sity of Dirac eigenvalues at zero, evaluated in the (large) 
volume V. Thus we expect that if the instantons and 
anti-instantons form a liquid, chiral symmetry will be 
broken, but that if they form instanton-anti-instanton 
molecules, chiral symmetry will be restored. 

If this interpretation is correct the chiral phase transi- 
tion should reflect itself also in a characteristic change of 
the localization properties of the lowest Dirac eigenvec- 
tors. This is the main motivation for our studies. 

In addition we are interested in comparing our lat- 
tice results with predictions from random matrix the- 
ory (RMT). Such a comparison allows us to distinguish 
generic from QCD-specific properties. Similar studies 
were performed before in Ref. In this paper, however, 
the authors analyzed a specific instanton-liquid approxi- 
mation to QCD while we study lattice QCD. 



As usual, we simulate a finite-temperature system by 
working on a lattice with Lt < Ls, where Lt (Ls) is the 
temporal (spatial) extent of the lattice. All boundary 
conditions are periodic except for the temporal bound- 
ary conditions of the fermions which are antiperiodic. 
The temperature is given by aT = 1/ Lt with the lattice 
spacing a. We fix Lt = 6 and vary a (and hence T) by 
changing (3 = Q/g^- 

We work in the quenched approximation with stag- 
gered fermions, i.e. with the Dirac operator 



D 
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where a^{x) = (_i)^i+ - +a^M-i and the are the link 
variables. The eigenvalues of D come in pairs of ±A, 
so we can restrict ourselves to positive eigenvalues in the 
following. In the continuum limit this action corresponds 
to four quark flavors. From now on we set a to 1. 

Quenched QCD has an additional symmetry of the 
gauge sector, which is spontaneously broken in the de- 
confined phase In the confined phase the expec- 
tation value of the Polyakov loop P is zero, whereas 
in the deconfined phase \P\ acquires an expectation 
value, and the phase of P clusters around the values 
6'p=arg(P)=0,±27r/3. The fermion action does not share 
the Z3 symmetry, so fcrmionic quantities can depend on 
the Z3 sector. It was found in Ref. that the chi- 
ral condensate computed from the configurations with 
Op — Q vanishes above Tc as expected. For Op = ±27r/3 
it remains finite in a certain temperature range above 
Tc- This behavior can be understood in Nambu-Jona- 
Lasinio models [||j7) and in RMT g. The point is that 
the boundary conditions of the Dirac operator are not 
invariant under transformations, and for Op — ±27r/3 
the new boundary conditions lead to a decrease of the 
Dirac eigenvalues ||]. According to the Banks-Casher re- 
lation, Eq. (|l|), this implies that the condensate will dis- 
appear only for larger temperatures, i.e. the transition 
temperature is higher for these sectors. For phenomeno- 
logical comparisons one should use the 9p=Q sector as 
this one is energetically favored for dynamical fermions. 
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(The sectors 6p = ±2tt/3 are physically equivalent to 
each other.) 

Large temperatures correspond to a small extension 
of the lattice in the temporal direction. The periodic 
boundary conditions for the gauge fields thus imply that 
the field equations have solutions with nontrivial topol- 
ogy which look rather like an instanton chain Such 
configurations, called calorons, have been the topic of in- 
tense investigations [Q. 

The index theorem tells us that the continuum Dirac 
operator has a chiral zero mode in a caloron field, so we 
will search for calorons by looking at the chiral and lo- 
calization properties of the eigenvectors of the low-lying 
eigenvalues of D. To calculate the eigenvalues and eigen- 
vectors we used the Arnoldi method as implemented in 
Ref. jl^ . This method allows us to choose the number of 
the lowest eigenstates and eigenvectors to be calculated. 

In the continuum a zero mode has a definite chiral- 
ity (= expectation value (75)) of ±1, while modes with 
A 7^ have (75) = 0. Staggered fermions possess only a 
restricted chiral symmetry, so in this case we expect (75) 
values between —1 and 1. In Fig. |l| we show a scatter 
plot of (75) against Dirac eigenvalue, for a temperature 
slightly above Tc. We can see that in all sectors the data 
points form clusters. There are some eigenmodes with 
very small eigenvalue, and (75) « ±0.2, while the bulk of 
the eigenvectors have a larger A, and (75) near zero. The 
modes labeled 1 to 3 in Figs. |l| and ^ will be discussed 
below. 

A large value of (75) suggests that the corresponding 
eigenvectors have a topological origin, in which case they 
should be associated with specific localized states. This 
interpretation is supported by the fact that these modes 
are approximately fourfold degenerate as they should be 
because of the flavor symmetry of staggered fermions in 
the continuum limit. One of our main results is that 
these states may be important for understanding the dif- 
ferences between the Z3 sectors. 

As a measure of the localization of our quark eigenvec- 
tors ipxix) (A is the Dirac eigenvalue, a a color index) we 
use a gauge- invariant inverse participation ratio (IPR), 
defined by 
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where V is the number of lattice sites and p\(x) is the 
gauge-invariant probability density 
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For a completely delocalized state (all p\{x) the same) 
one finds /2 = I, whereas a state localized on a single 
lattice site (only one non-zero p\{x)) would have I2 = 
V. (The staggered fermions' chiral symmetry means that 
in fact an eigenstate can never be completely localized, 
since all eigenstates must have half their probability on 
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FIG. 1. (75) vs. A. The measurements are on a 12^ x 6 
lattice at /? = 6.10 for both sectors, 6p = and 9p — ±2n/3. 
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FIG. 2. A scatter plot of the localization vs. chirality. 

even sites, and half on odd.) I2 is also a measure of the 
standard deviation of px (x) , 



/2 - 1 = 



VpI 



(5) 



Chiral RMT predicts that the average value of I2 is 
[N, + \)V V 



{l2) = 



Ac\^ + 2 



I + -L forA^, >3. (6) 



To elucidate the relevance of the IPR let us note that in 
condensed matter physics the size of the IPR decides e.g. 
whether a disordered mesoscopic sample is a metal {I2 
close to the RMT prediction) or an insulator {I2 large). 

We have investigated the behavior of I2 on finite tem- 
perature lattices, on both sides of the deconfinement 
phase transition, which lies at (3 ~ 5.89 for it = 6 in 
the thermodynamic limit |ll|. 

Slightly above Tc the localization properties show char- 
acteristic differences for the different sectors. In 
the real sector the effect of localization is strongly pro- 
nounced while it is nearly absent in the 9p = ±27r/3 
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FIG. 3. From top to bottom, the average IPR for temper- 
atures below, slightly above, and well above the chiral phase 
transition plotted separately for the Z3 sector with 9p = 
and the Z3 sectors with 9p = ±27r/3. The average was per- 
formed for the 2°"*, . . . , 20*'' eigenvalue of each gauge field 
configuration. The dashed line is the RMT prediction, 4/3. 

sectors where chiral symmetry is still broken for these 
temperatures. This is illustrated in Figs. |^ and ^ A dif- 
ference between localization in the different sectors tells 
us immediately that the eigenstates must be extended 
in the time direction, because an eigenstate localized in 
time would be unaffected by a Z3 transformation. 

We studied these characteristic differences for high 
temperatures in more detail and found features which 
suggest that they can be attributed to the effects of 
calorons, as we shall now explain. 

Harrington and Shepard found exact SU(2) instan- 
ton solutions at finite temperature M, which they called 



calorons (for a review see Ref. [Q). Their solution con- 
sists of a one-dimensional chain of instantons, with a pe- 
riodic repeat at a distance 1/T. There is a single physi- 
cally relevant parameter, RT, which gives the instanton 
radius R, relative to the scale set by the temperature. 
When RT ^ 1 the solution looks like an isolated instan- 
ton, but when RT 1 finite temperature effects become 
significant. This instanton chain satisfies the 't Hooft 
ansatz so the techniques described in Ref. can 
be used to find the fermionic zero modes |l5|] . For small 
RT these look like the zero modes of an isolated instan- 
ton, but as RT increases, the modes become extended 
in time, though they remain localized in space. These 
SU(2) solutions can easily be embedded in SU{3), lead- 
ing to solutions which correspond to the 9p = sec- 
tor. To investigate the other sectors, we need solutions 
where the time-like Polyakov loop has a non-zero phase at 
space-like infinity. We can construct such solutions either 
by adding a constant A4 background field with a color 
that commutes with the SU(2) subgroup containing the 
caloron, or by choosing the fermion field boundary con- 
dition appropriately. In either case the fermionic zero 
modes are readily constructed. One can also find pure 
SU(2) solutions with a Polyakov loop background [ p^ . 
It would be interesting to study the relationship between 
these solutions and the solutions found by embedding the 
Harrington-Shepard solution in SU(3). 

The localization of the embedded caloron's fermion 
zero mode depends heavily on the Z3 sector. Asymp- 
totically the modes fall off like 

IVI' ocexp[-2(7r- |0p|)rr]/r2 (7) 

where r is the three-dimensional distance from the 
caloron axis. (At 9p — 0, this agrees with the behav- 
ior given in Ref. (l2j.) This is unlike the case of a single 
instanton, which has \^p\'^ which drops off like a power of 
r. We see from Eq. (|^) that the correlation length for 
the zero mode is smallest in the real sector {dp = 0), and 
that the complex Z3 sectors {Op = ±27r/3) have modes in 
which the radius is about three times greater, and so they 
occupy a much larger volume. It is thus tempting to as- 
sume that the strong difference between the localization 
properties in the Z3 sectors is due to the fact that our 
lattice configurations contain calorons. If so, the local- 
ization in four dimensions should show the characteristic 
string-like pattern discussed above. 

In Fig. Q the localization itself and the local density 
of the expectation value of 75 are shown for a highly lo- 
calized state in the 9p = Q sector (mode 1 of Figs. |l| 
and H). To plot the complete four-dimensional lattice 
we have introduced the two coordinates i — x + 12t 
and j — y + 12z with the Euclidean lattice coordinates 
X, ?/, z = 0, 1, 2, . . . , 11 and t — 0, ... ,5. This coordinate 
system represents the lattice as a 6 x 12 array of x-y 
slices, which generates the approximately periodic struc- 
ture visible in the plot. From the upper plot, we see that 
the state is indeed spatially localized, but extended in 
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FIG. 4. A 'caloron' state (mode 1 of Figs. ^ and ^) on a 
12^ X 6 lattice for jS = 6.1 in the 6p ^ sector. We have 
introduced the coordinates i — x + 12t and j = y + 12z with 
the lattice coordinates x,y, z = 0, 1, . . . , 11 and t — 0, ... ,5. 
The gauge- invariant density is plotted above, the expectation 
value of 75 below. Both show localization in space but not in 
time. 



time. The lower plot demonstrates that the eigenstate 
in question is an approximate chiral eigenstate. A con- 
tinuum caloron would possess an exact chiral eigenstate. 
Unfortunately, staggered fermions exhibit only part of 
the continuum chiral symmetry, which is the reason why 
we do not find perfect 75 eigenstates. It would therefore 
be very interesting to repeat these studies with Ginsparg- 
Wilson fermions which have much better chiral proper- 
ties. 

The corresponding pictures for mode 2 look similar to 
Fig. ^. The width of the mode is, however, noticeably 
larger than that of mode 1. Mode 3 serves as an exam- 
ple for a mode with low (75), high I2, and an eigenvalue 
which lies in the bulk of the spectrum. Although (75) is 
small, the (75)-density has a region of large positive val- 
ues right next to a region where the density is large and 
negative. This suggests mode 3 might be an instanton- 
anti-instanton molecule. 



Modes related to instantons and calorons were rece ntly 
investigated in lattice QCD with dynamical fermions . 
While the results below Tc are consistent with expecta- 
tions based on an instanton liquid picture, the interpre- 
tation of the observations above Tc appears to be less 
straightforward. 

Let us summarize: We have analyzed the localization 
properties of quark eigenstates in quenched lattice QCD. 
By concentrating on the low eigenvalues we could char- 
acterize semiclassical properties of the gauge field config- 
urations without any cooling. (For investigations using 
cooling, see e.g. Ref. 10.) For temperatures above Tc 
we found isolated modes with definite handedness which 
show all the properties of fermion states associated with 
calorons, in particular localization in space but not in 
time. 

We wish to thank A. V. Belitsky, Ph. de Forcrand, 
and T. Schafer for useful conversations. 
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